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Abstract 

In this paper we derive boundary integral identities for the bi-Laplacian 
eigenvalue problems under Dirichlet, Navier and simply-supported bound- 
ary conditions. By using these integral identities, we first obtain unique- 
ness criteria for the solutions of the bi-Laplacian eigenvalue problems, and 
then prove that each eigenvalue of the problem with simply-supported 
boundary conditions increases strictly with Poisson's ratio, thereby show- 
ing that each natural frequency of a simply-supported plate increases 
strictly with Poisson's ratio. In addition, we obtain boundary integral 
representations for the strain energies of the vibrating plates under the 
three boundary conditions. 
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1 Introduction 

The bi-Laplacian boundary value problems with Dirichlet, Navier and simply- 
supported boundary conditions are classical boundary value problems. Among 
these problems, Dirichlet and simply-supported problems are of more impor- 
tance in solid mechanics. In this paper we shall first study the eigenvalues 
of the three problems, and then pay special attention to the simply-supported 
problem, in contrast to many publications which only focused on the other two 
problems. In fact, the Navier boundary conditions are a special case of the 
simply-supported boundary conditions and the Dirichlet (clamped-edge) prob- 
lem can be approximated by the simply-supported boundary value problem if 
the interior domain is far enough from the boundary. 

It is known that resonance problems are a major concern in mechanics. Res- 
onant frequency of a solid structure is an instrumental parameter that greatly 
affects the structural dynamic behaviors governed by the solid mechanics the- 
ory [1], [!]■ In this section, we shall simplify the simply-supported boundary 
conditions and give the problem statements for the three types of the eigenvalue 
problems. F. Rellich introduced a trial function X]fc=i ^k-§^ = x ■ \7u, applied 
it to the Dirichlet Laplacian operator and obtained a boundary integral identity 
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for the Dirichlet Laplacian problem in 1940 [3]. This idea has been generahzed 
and appHed to elliptic PDE problems [4j, [5], [6]. In order to show the depen- 
dence of the eigenvalues on the parameter in the boundary conditions of the 
vibrating plates, we shall follow Rellich's idea to derive three boundary integral 
expressions for the Dirichlet, Navier and simply-supported bi-Laplacian prob- 
lems and shall establish uniqueness criteria for the solutions of the eigenvalue 
problems in Section 2. Based on these integral identities for the bi-Laplacian 
problems we show how the eigenvalues of the simply-supported problem and 
the resonant frequencies of the vibrating plate are influenced by Poisson's ratio 
through the new uniqueness theorems in Section 3. In addition, we derive three 
boundary integral expressions for the strain energies of the vibrating plates at 
resonance using the boundary integral identities. These boundary integral ex- 
pressions can be used to calculate the total strain energies when only boundary 
value data of the plates are available. 

1.1 Dirichlet and Navior Bi-Laplacian Problems 

Assume that 57 C > 2, is a bounded domain having a C^'^ boundary 

9f2 (0 < /3 < 1). Let A be an eigenvalue for which the Dirichlet eigenvalue 
problem, 

AAU ^ AU in f], 

u = |v;7| on an, (i) 

has a nontrivial solution, and let A be an eigenvalue for which the Navier eigen- 
value problem, 

AAV" = XV in Q, 

V = A\/ = on on, (2) 

has a nontrivial solution, where AA = A^ denotes the n-dimensional bi-Laplacian 
and U and V are respectively the eigenf unctions of Problems H]) and 

1.2 Simply supported Boundary Conditions 

The supported boundary conditions (short for simply supported boundary con- 
ditions) in their standard form obtained from their mechanical property and 
implementation manner have been applied to the bi-Laplacian problems [I] , [5] . 
These conditions are natural, based on the physical ground. In this subsection, 
we shall first simplify the standard form, and then define the corresponding 
eigenvalue problems. 

We consider a thin elastic plate which is assumed to be homogeneous and 
isotropic under the supported boundary conditions: 

Wlen = 0, 
M^lon = 0, 
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where W is the deflection of the plate in the vertical direction and is the 
bending moment with respect to the v direction. These boundary conditions in 
two dimensions can be expressed as (See [TU], Expression (110), p. 94.) 



M^lan = ^AW^Ian + (1 - m) 



COS 9——: — h 2sin9cos9— — - — f- sm 



dx^ dxdy dy^ 



where /x is Poisson's ratio, an elasticity constant (0 < /i < 1), the normal unit 
vector V is directed outward from f2 and 9 is the angle between the normal v 
and X-axis {cos9 = cos{x,v) = v^)- It is shown that the factor in the square 
brackets of the above equality is equal to (See the appendix to this 

paper). This can be written as 

cos + 2sm9cos9——- — h sm 9 = , on oil. 
OX'' oxoy oy' ov' 

Thus, we obtain 

M^\dn = ^il^W\^n + (1 - i^)-Q^\oa = ^■ 

Using the above definition and the relation AVt^lgji — + (n — l)K-^-|an for 
n — 2, we obtain 

dW dW 
AW\an - (1 - = g^lon' 

where k is the mean curvature of dfl and cq = (1 — /i)K > for convex domains. 

1.3 The Supported Bi-Laplacian Problem and Basic Iden- 
tities 

We are now in a position to define the supported bi-Laplacian eigenvalue prob- 
lem. Let J7 C i?^ be defined as in Subsection 1.1 and let 7 be an eigenvalue 
for which the bi-Laplacian eigenvalue problem under the supported boundary 
conditions, 

AAT4^ = -/W in Q, 

w ^ on on, 

dW 

AW = Co— on dVL, (3) 
ov 

has a nontrivial solution. Elliptic regularity theorems ensure that Problems ([1]) , 
([2]) and ([3]) will respectively have nontrivial solution U , V and W G C^{Vt). 
Assuming that a and u are the eigenvalue and the corresponding eigenfunction 
of Problem ([T]) or ([2|) or ([3]), Green's identities show that 

I {au)udn = [ {AAu)udn = [ ^^^udS - [ VAu ■ VudD, 
Ja Jn Jan Jq 

^ [ ^udS~ I ^AudS+ f {Aufdn. (4) 
Jan oi' Jan ov 
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The common boundary condition of Problems (P), ^ and ([3]), U\dn = V\gn = 
W\an = 0, leads to f^l,^ = |VC/||,^, |^|,^ = |Vy||g^ and ^\^^ = 
|VW^||gj^. For Problems H]) and ([2]), both the first and second surface inte- 
grals on dft vanish. Therefore, the eigenvalues of the two problems can be 
expressed by the same formula 

^^'^ = r In ^ 
Jo ^ dn 

where u can be U or V. For Problem ([3]), the identities shown in (U) shows that 



(6) 



The positiveness of the eigenvalues are proved by the theory of elliptic partial 
difi^erential equations [6], [7], [8]. 



2 Boundary Integral Expressions for 
Bi-Laplacian Problems 

Theorem 2.1 Let 57 be defined as in section 1.1 and let J7 be a nontrivial 
solution of Problems H]). Then, Let fl be defined as in section 1.1 and let U be 
a nontrivial solution of Problems ([1]). Then, 

Theorem 2.3 Let 51 be defined as in subsection 1.2 and let ^ be a nontrivial 
solution of Problem ([2]). Then, 

Theorem 2.5 Let ft be defined as in subsection 1.3 and let be a nontrivial 
solution of Problem Then, 

, _ - Uni- ■ -)ico^?dS + 2 J,^ix ■ ^m^rdS 
Proofs of theorems 2.1, 2.3 and 2.5 

We first multiply both sides of the bi-Laplacian eigen equation of Problems 
([!]), (12) and ^ by the test function {x ■ Vu), and integrate over 51 to obtain 

/ {AAu){x ■ Wu)dn = / au{x-Vu)dn. (10) 
Jn Jn 
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Applying Green's theorem to piT)) . the right-hand side of pUj) becomes 
Jn dxk 2 ^ dxk 2 2 Jgf^ 



fc=i 

and the left-hand side of ITUl) becomes 



/ (AAw)(2:-Vu)df7 = - / \7Au-y{x-\^u)dn+ [ ^^{x-yu)dS. (12) 

The (generalized) volume integral term of the right-hand side of can be 
written as 

VAu • V(x • Vu)ai2 = - / > T, —dil 

dAu I du d'^u \ 

^ dxi \ dxi ^ ^ dxidxj ' 




where Vj — cos{xj^v). The volume integral term of the right-hand side of the 
above equality can be further written as 




From (HI, we have 

a [ u^dn= [ ^^udS- [ ^AudS+ [ {Aufdfl. 
Jn Jan ov J an dv 

Expressing J^{Au)^dil in terms of the integrals in the above equality, (|13p 
becomes 

{2-l)(al u^dil- I ^udS+ I ^Auds]+l [ {Auf{x-y)dS 
2 V 7a Jon dv Jg^ dv ) 2 Jg^ 

(14) 
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Substituting into and equating the result to pUj) . the following 
identity results 

(2 - ^)(a[u^dn-f ^udS+ f ?^Auds]+l[ {x ■ ,,){Au)^dS 
2 V Jon '^^ J do. '^^ / 2 Jq^ 

. I du ^ d'^u \ ,„ /" dAu 
Au\— + } xiiyj- — ^ \ dS + — — (x-Vu)dS 
on ^ dxidxj J Jq<^ dv 

= /■ + ^ / {x-v)u'dS, 

where the relation 

/ rP-ii \ t I rill rP'ii \ 

dS 



O u 

'3 



is used. After cancelation, the above identity can be written as 

2a / u'^dn -- / {x ■ i^){AufdS + / Au — +Va;jzy 

Jan 9iy 2 Jg^ 2 Jq 



^ dxidxj 



dAu 9u\ ,^ 

— — u - Au— I db 
ov ou 



We first consider the Dirichlet boundary conditions of Problem Let u = U, 
a = A and U\ga — \ VU\\gQ — 0; then ([TS]) becomes 

2a f U^dn ^ -I f {AUfix-v)dS+ I AU x^v^^^^^\ dS. 
Jn ^ Jan Jon ' dx.dx, j 

Solving for A using ^'^ ^ = i^ii^iTr-^ on dfl, then 



A 



hni^ ■ ^){^UfdS + 2 AC/ (Er, aSm^) dS 



^ -'Iani---)&'dS + 2f^,g^{ix..)^) 

Aj^u^dn 

Therefore, 



dS 



A = 



Aj^u^dn 



Based on the fact that the eigenvalues of Problems (1) are greater than zero, 
we obtain 



Corollary 2.2 Problem (1) admits only a trivial solution U £ C^(f2) for 
C i?",n > 2 if 

= ^lao = and = . 

Secondly, we impose the Navier boundary conditions of Problem ([2]). Let u = 
a = X and V\dn = AV\dn — 0; then (fTS]) becomes 



Jan Jon 



I an Jon 
Solving for A, we obtain 

Based on the fact that the eigenvalues of Problems (2) is greater than zero, we 
obtain 

Corollary 2.4 Problem (2) admits only a trivial solution V G C'^{Cl) for 
n C > 2 if 

fdAV dV\ 
V\9n = AV\an = and or {x ■ v)— j j^^ = 0. 

Since AV\an = + (n — l)K^\Qn = 0, the if conditions in Corollary 2.4 
can be written as 

d'^V , dV . 

V\an = and = -{n - and 

/d^V d^V . dV . 

-TT^ + (n — -— - =0 or (x ■ v) — — = 

Next, we impose the supported boundary conditions of Problem (|3]). We 
assume the domain J7 is convex, which is the common case for plates. Letting 
u = W,a = ^,n — 2 and inserting W — 0|ao and AW — co^|an into (|T5|l . 
we obtain 

27 / VF^df] = / Co— 2^5+ / cq— Va^.^^jTr^ 



/ ^U'-^)(^-^w)ds, 

Jan 9v \ ov J 
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with Co > 0. Solving for 7, we obtain 



where the relations g^.g^ . = t^i^^j^j^ and a; • VW = ^-{x ■ v) on d^l are used. 
This completes the proof of Theorem 2.5. 

Based on the fact that the eigenvalues of Problems (3) are greater than zero, 
we obtain 

Corollary 2.6 Problem (3) admits only a trivial solution W G C'^{Vl) for 
C i?2 if 

Corollary 2.7 Corollaries 2.2, 2.4, 2.6 illustrate that a nontrivial eigenfunc- 
tion cannot simultaneously be an eigenfunction of the Dirichlet and supported 
(or Navier) plate of the same shape. 

3 Application to Solid Mechanics 

Plates (short for elastic uniform-thin plates) are a particular two-dimensional 
representation of a three-dimensional solid, which have a much smaller thickness 
in comparison with the in-plane dimensions [1], [2]. In the following, we will 
apply our theorems and corollaries to vibrating plate problems. 

3.1 Effect of Poisson's Ratio on Eigenvalues 
of Supported Problem 

Theorem 3.1 Let 17 be a convex domain defined as in section 1.3 with k{x) > 
0, Va; C dfl and let W{x,Iji) € Ci'^{Cl, (0, 1)) be a nontrivial solution of the 
supported eigenvalue problem. Then, 

where W{x,fi) e C^''^ {Cl, (0,1)) denotes W{x,-) e C^'>^{n) and W{-,fi) G 
C\0,1). 

Proof of Theorem 3.1 By using AAW — jW of Problem ^ and differ- 
entiating the identity with respect to n, we have ^^^^ — 7^ + T^f^- 
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Since is continuously differentiable, we obtain 



Therefore, 



dAAW _ _^_^dW _ dW ^dj 



dW dW n 9t 

WAA'^ = ^W^ + W^^ (17) 
Ofi dfi Ofi 



Integrating both sides of the above identity in Q. and using the Green's identities, 
the left-hand side of (fT7|) becomes 

dW r d dW r dW dW r d'^W 

WAA^dn ^ / W^A^dS- / ^A^d5+ / ^^--^AWdS 
OjJ. Jon ou dfi Jq^ du dfi Jg^^ dfidv 

dAW dW f ^ ^„,dW ,^ 
-dS+ / AAW-—dn 



I an dv dfi dfi 

Using AAW^ = 7W and identity p7|) . and imposing the supported boundary 

dw I _ „ aw_\ 
dv Ian ^0 dv Ian' 



conditions M^|an = and AW\Qn = (1 - m)k7^L„ = cotj^L^, we obtain 



r.r..dW,^ f dW dW r dW d^W ,^ 

WAA—d^l = 1-^ +co—^)dS+ Co ^ dS 

dfi Jan ov \ dv dvdfij Jg^^ dv dvdfi 

f dW f ,dW,2 f dW 

+ / iW^dS= ^^fdS+ Hw^dQ, (18) 
Jan dfi Jgn 'dfi' Ja dfi 

wliere = and ^\g^ = (k^ + co^)\g^ is used. By using the 

identity WAAW^dQ = ^ (w^^ + 7^^^) dn, ^ can be written as 



Ids 



k\^—\ dS = / W'^-r-dQ 



dv Jq dfi 



/an 
Hence, 

d-f L^Kl^f dS , , 

dfi- J^w^dn ^^'^ 

which shows that 7 is a nondecreasing function of fi. 

Furthermore, for those domains with k > on dCl, from (jl9p . we obtain 
that if |l = 0, then ^ = on 917. Applying Corollary 2.6, if C i?^ with 
k{x) > on dil, then = in fi. This contradicts the assumption that W is 
a nontrivial solution of Problem ([3]). Hence, |j ^ 0. Therefore, 

dfi 

for all convex domains with k > on dfi. This completes the proof of theorem 
3.1. 
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Corollary 3.2 Each eigenvalue of simply-supported convex plates with k > 
is strictly monotonic in fi. Equivalently, a nontrivial eigenfunction cannot 
simultaneously be an eigenfunciton of two simply-supported plates of the same 
shape with different values of jj, if the boundary of the plate is convex and there 
is no straight segment on the boundary (k > 0). 

It is known that a natural (intrinsic) frequency of a vibrating plate uj = 

7^, where m denotes the equivalent mass per unit area of the plate, D = 

j2 {i-ii'^) ' ^ denotes the modulus of elasticity and h is the thickness of the plate. 
Therefore, Theorem 3.1 leads to the following corollary. 

Corollary 3.3 Each natural frequency of simply-supported convex plates 
with K > increases strictly with Poisson's ratio. 

The property of the internal parts of the plates far enough from the boundary 
can be approximated by the clamped boundary conditions, and thus the result 
is applicable to the Dirichlet problem in the internal parts. 



3.2 Boundary Integral Identities for the Eigen-Problems 

The Strain energy of a mechanical structure provides a good measure for ex- 
ceeded stresses and strains or exceeded strain energy in the structure. It is often 
used, with certain failure criteria of materials, to evaluate if a structure under 
certain loads is in a safe condition. Identities representing the strain energy 
conservation of elastostatic problems have been developed and used to prove 
the uniqueness of the equilibrium solutions in elastostatics [12] , [H] . 

For the Dirichlet Problem ([1]), applying Theorem 2.1 and (O, we have 

U^u?dn /^^(...)(0)^d5 

Multiplying both sides of the above identity by U^dfl, the following identity 
yields 

/ (AUfdn^l I {x-v){^fdS (20) 
For the Navier Problem ([2]), applying Theorem 2.2 and (O, we have 

J^iAvfdn _ ~Sg^{x.vv)'-^ds 

Multiplying both sides of the above identity by V^dfl, the following iden- 
tity yields 

(AVfdn^-l f (x-VV)^^dS (21) 
2 Jdn ov 

For Problem ([3]), applying Theorem 3.5 and ©, we have 

UAWfd^ - co(^)^d5 - !,^{x ■ u){c,^-^fdS + 2jg^ix ■ ^)%^i^)^dS 
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where co{x) = k(.t)(1 — fi). Multiplying both sides of the above identity by 
W^dH., the following identity results 



an I 



on dv 



+ 



Jan oi' 



(22) 



3.3 Boundary Integral Expressions for Stain Energies 

In the following we shall derive the boundary-integral expressions for the strain 
energies of the vibrating thin plates. 

The general expression for the strain energy of a bent thin plate is given as 
(See Hn], p.95.) 



-D 



dx^ dy"^ dxdy 



dn, (23) 



where u is the deflection of the plate in the vertical direction and the integration 
is extended over the entire surface of the plate. Here u can be the solutions of 
Problems H]), ([2]) and ([3]). In the derivation of the strain energy expressions the 
following identity is needed (See [13], p. 87) 



dx^ dy^ dxdy 



dQ 



an 



du d"^ 



1 ,du 



1 ,du^ 



dS, 



with K > 0. Under the Dirichlet boundary conditions this identity reduces to 



-? 



dfl = 0, 



(24) 



dx'^ dy^ dxdy' 

Under the Navier or supported boundary conditions this identity reduces to 

dQ = l [ -i^fdS. (25) 

For the Dirichlet Problem, setting n — 2 and substituting (PD|) and into 
we obtain 



d'^u d'^u d^u 
dx"^ dy"^ dxdy 



1 

-D 



an 



(26) 



where is the strain energy of the resonant plate under the Dirichlet boundary 
conditions. 

For the Navier Problem, setting n = 2 and substituting (PT|) and (pS)) into 
(l23l) we obtain 



E: 



N 



D 



an 



1, .dVdAV ^ 1,9^,2' 

2 dv du K dv 



dS 



(27) 
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where is the strain energy of the resonant plate under the Navier boundary 
conditions. For the supported plate problem, setting n — 2 and substituting 
(1221) and (12S1) into ^ we obtain 

8 Jq^ dv 4 Jg^ dv dv 

Z? /• ,dW . D f l,dW , 

+ T/ (^o{-7T-) dS ~ —{1 - n) -{ — ) dS, 
2 Jdn 2 Jq^ k dv 

where Ef is the strain energy of the resonant plate under the supported bound- 
ary conditions. The supported problem can be modeled as a hinged plate at- 
tached to the boundary with a torsion spring. When the spring constant goes 
to infinite, the problem becomes a clamped boundary plate. 

Remark 4.1 The above three forms for E^ , E^ and E^ are the boundary 
integral expressions for the strain energies, which simplify the calculation of the 
strain energy in a plane area to the boundary of the area at resonance. These 
expressions can be used to simplify the calculation of the strain energies if only 
boundary data are available. 
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4 Appendix 1 

We shall prove the equality, cos^O^ + 2sin6cos6^ + sin^9^ = ^ on 
951, used in simplifying the supported boundary conditions. 

Let (si, . . . , Sn-i) be a local coordinates on dQ, i.e., dV, is locally represented 

by 

X : K"-^ -> dn 

Let 

v{si, . . . , Sn-l) 

be its unit outer normal of 951 at X{si, . . . , s„_i). We also define a local coor- 
dinate system near dfl by 

X{si, . . . , S„-l, i) = X{si, . . . , Sn-l) - tv{si, . . . , S„_l). 

In two dimensional case, for any u defined near 951, we choose the coordinate 
system (s, i), such that t is the distance from a point x in 51 to dVl and s is the 
arc-length parameter on 951 in the counter-clockwise direction. On 951 we write 
V = (vx, Vy) and have 

Vt = -V and Vs-Vt^ 0. 
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Hence on dfl. 



d'^u d'^u ( dt ds dt ds\ d'^u dt dt du d'^t 



dx^ dtds \ dx dx dx dx J dt^ dx dx dt dx"^ 

Let cos9 = cos{x, — v^, sinO — cos{y, v) = Vy and define 

= — -^cos y + 2 cosOsmO + -—^sm 
ox-' oxoy oy^ 

Substituting ([25)1 into the above expression for G, then 

where we have used the fact that 

d'^t d f dt dx dt dy 



dt'^ dt \dx dt ' dy dt 



dH dH 2 dH 

—r + 2V^U 

dx'^ 

This complete the proof of the equality. 



^^'^d^^^'^^^yd^y^^yW^-""' 



References 

[1] Timoshenko, S., D.H. Young and W. Weaver, JR. Vibration Problems in 
Engineering, John Wiler & Sons, 1974. 

[2] Timoshenko, S., Theory of Plates and Shells, (second edition). Maple Press 
Company, York, PA, 1979. 

[3] Rellich, F., Darstellung der Eigenwerte Au + Am durch ein Randintegral, 
Math. Zeit.,Vol. 46, 1940, pp. 635-646. 

[4] Payne, L.E., Inequalities for Eigenvalues of Supported and Free Plates,^. 
Rational Mech. Anal. 4,(1955), pp. 517-529. 

[5] Payne, L.E., Inequalities for Eigenvalues of Membrane and Plates, Quart. 
Appl. Math. 16, 1958, pp. 111-120. 

[6] Drabek, P., Global Bifurcation Result for the p-biharmonic Operator, Elec- 
tron. J. Differential Equations, 2008, No.48, pp. 19-23. 



13 



[7] WangjM., Fu, A. and Han, Y., Lower Bounds for Dirichlet Eigenvalues 
of Higher-order Elliptic Operators, MATHEMATICA APPLICATA,2000, 
13(4), pp. 21-24. 

[8] Payne, L.E. and Weinberger, H. F., Lower Bounds for Vibration Frequencies 
of Elastically Supported Membranes and Plates, J. Soc. Indus. Appl. Math. 
Vol. 5, No. 4, December, 1955, pp. 171-182. 

[9] Ranch, J., Partial Differential Equations, Springer- Verlag, New York, 1991. 

[10] Timoshcnko, S., Theory of Plates and Shells, (First Edition), Maple Press 

Company, York, PA, 1940. 

[11] Mazya, V.G. and Nazarov, S. A., Paradoxes of Limit Passage in Solutions 
of Boundary Value Problems in Solving the Approximation of Smooth Do- 
mains By Polygonal Domains, Math. USSR Izvestiya, Vol. 29, No. 3, 1987, 
pp. 511-533. 

[12] Knops, R.J., Uniqueness and Continuous Data Dependence in the Elastic 
Cylinder, International Journal of Non-Linear Mechanics, 36, 2001, pp. 489- 
499. 

[13] Birman, M. sh. Variational Methods for Solving Boundary Value Problems 
Analogous to Trefftz' Method, Leningrad University, A. A. Zhdanova, 1956, 
pp. 69-89. 

[14] Knops, R.J. and Stuart, C.A. Quasiconvexity and Uniqueness of Equilib- 
rium Solutions in Nonlinear Elasticity, Arch. Rational Mech. Anal. 86/3, 
1984, pp233-249. 

[15] Friedlander, L Remarks on the Membrane and, Buckling Eigenvalues for 
Planar Domains, Moscow Mathematical Journal, vol. 4, no. 2, April- June, 
2004, pp. 369-375. 

[16] Robert, van der Vorst, Fourth order elliptic equations with critical growth, 
C.R. Acad. Sci. Paris, t. 320, Serie I, 1995, pp. 295-299. 



14 



